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We explore spin- 1/2 triangular antiferromagnets with both easy-plane and lattice exchange anisotropics by 
employing a dual vortex mapping followed by a fermionization of the vortices. Over a broad range of exchange 
anisotropy, this approach leads naturally to a "critical" spin liquid — the algebraic vortex liquid — which appears 
to be distinct from other known spin liquids. We present a detailed characterization of this state, which is de- 
scribed in terms of non-compact QED3 with an emergent SU(4) symmetry. Descendant phases of the algebraic 
vortex liquid are also explored, which include the Kalmeyer-Laughlin spin liquid, a variety of magnetically or- 
dered states such as the well known coplanar spiral state, and supersolids. In the range of exchange anisotropy 
where the "square lattice" Neel ground state arises, we demonstrate that anomalous "roton" minima in the 
excitation spectrum recently reported in series expansions can be accounted for within our approach. 

PACS numbers: 75.10.Jm, 75.4().Gb 



I. INTRODUCTION 

A fundamental theoretical challenge in strongly correlated 
systems lies in understanding the behavior of frustrated quan- 
tum magnets, whose properties often bear little resemblance 
to those of their classical counterparts. In the most exotic 
scenario, quantum fluctuations are sufficiently strong to dis- 
order the system even at zero temperature, and a spin liquid 
ground state emerges. Historically, Anderson originally sug- 
gested that the spin- 1/2 Heisenberg triangular antiferromag- 
net may realize such a quantum-disordered ground stateii It 
is now recognized that with only nearest-neighbor exchange 
the true ground state on the triangular lattice is the magneti- 
cally ordered \/3 x \/3 phase, though the order is significantly 
diminished relative to the classical state. ^ ■^ It is conceivable, 
then, that a spin liquid may arise with not too drastic perturba- 
tions to the model, and the triangular lattice has thus remained 
a prominent setting in the search for two-dimensional spin liq- 
uids. 

Recent experiments on the spin- 1/2 anisotropic triangu- 
lar antifeiTomagnet CS2CUCI4 stimulated renewed interest in 
possible spin-liquid phases proximate to the nearest-neighbor 
Heisenberg model. "^-^ This material is accurately modeled 
by an anisotropic Heisenberg Hamiltonian supplemented by 
a weak Dzyaloshinskii-Moriya interactioni^ Although long- 
range spiral order develops at temperatures T < 0.62K, the 
dynamical structure factor measured via neutron scattering ex- 
hibits "critical" power laws at intermediate energies, both in 
the ordered phase and in a range of temperatures above T^r. 
This unusual power law behavior in the excitation spectrum is 
highly suggestive of spinon deconfinement that is characteris- 
tic of spin liquids. 

A variety of theoretical approaches have been employed to 
capture spin liquids on the triangular lattice. Kalmeyer and 
Laughlin exploited a mapping between the spin- 1/2 Heisen- 
berg model and hard-core bosons in a magnetic field to obtain 
a "chiral" spin liquid which breaks time-reversal symmetry 
Their arguments were subsequently reformulated by Yang et 
al.^ who arrived at the chiral spin liquid by fermionizing 
the spins using Chern-Simons flux attachment'" and expand- 



ing around a "flux-smeared" mean-field state. Using a slave- 
boson representation of the spin operators, Sachdev explored 
an Sp{N) generalization of the Heisenberg model, and in the 
large-A^ limit obtained a Z2 spin liquid ground state, which 
breaks no symmetriesiii The Z2 spin liquid was later realized 
microscopically in a quantum dimer model on the triangular 
lattice.'^ Finally, a large class of spin liquids was studied by 
Zhou and Wen using a slave fermion representation of the 
spins. Whereas excitations in both the chiral and Z2 spin 
liquids are gapped, the slave fermion mean-field approach can 
give rise to so called "algebraic spin liquids", which admit 
gapless spin excitations and power-law spin correlations. 

In this paper we pursue an alternate approach to the spin- 
1/2 triangular antifeiTomagnet, and use vortex duality to at- 
tack the problem coming from the easy-plane regime. Duality 
has been a powerful tool for exploring unconventional phases 
such as valence bond solids and spin liquids in quantum spin 
systems, '''•'^•'^•'^•'** as well as complex charge-ordered states 
in bosonic systemsii2i22iSii22i22i2i The main difficulty here is 
that vortices are at finite density, which is familiar from dual 
approaches to the fractional quantum Hall problem. As an 
initial step towards applying duality to frustrated spin sys- 
tems, in Ref. 25 we examined integer-spin triangular antifeiTO- 
magnets with easy-plane symmetry from the vortex perspec- 
tive. By fermionizing the vortices using Chern-Simons flux 
attachment, it was shown that an effective low-energy dual 
formulation can be derived, which was argued to reproduce 
the physics of a more direct Landau-Ginzburg-Wilson analy- 
sis of the spin model. While this approach is reminiscent of 
the spin fermionization adopted by Yang et al.^ we empha- 
size that alternatively fermionizing vortices is advantageous 
because the vortices interact logarithmically, which allows for 
a more controllable treatment of Chern-Simons gauge fluctu- 
ations. 

Here we extend the fermionized vortex approach to the 
spin- 1/2 triangular antiferromagnet with easy-plane symme- 
try and anisotropic nearest-neighbor exchanges J and J' as 
shown in Fig. ^ This formalism allows us to explore the 
phase diagram of the spin model in a setting where a more 
conventional Landau analysis of the spin model is not acces- 



sible due to Berry phases. Remarkably, over a broad range 
of anisotropy J'/ J < 1.4 this approach leads naturally, with 
the simplest flux-smeared mean-field starting point, to a novel 
"critical" spin liquid that we will refer to as the algebraic vor- 
tex liquid. This state was introduced earlier and applied to 
CS2CUCI4 in a short letter, Ref. and is characterized in 
detail here. Schematically, vortices form a critical state with 
four Dirac nodes, and interact via a fluctuating gauge field 
representing the original boson current fluctuations. Already 
on the mean-field level, the gapless character of the vortex 
state implies power-law and spin correlations at spe- 
cific wave vectors shown in Figs.|6land0 respectively. Such 
momenta for low-energy excitations in the spin system are de- 
termined by short-distance physics in the frustrated magnet, 
and we propose that this physics is well-captured in the vortex 
treatment. 

Going beyond a mean-field analysis, we argue that the al- 
gebraic vortex liquid is described by QED3 with an emer- 
gent global SU(4) flavor symmetry, which has further im- 
plications for the dynamical spin correlations. In particular, 
as a consequence of the SU(4) symmetry the in-plane spin 
structure factor exhibits enhanced universal power law cor- 
relations with the same exponent at several momenta in the 
Brillouin zone: the spiral ordering wave vectors ±Q and mo- 
menta Ki 2,3 at the midpoints of the Brillouin zone edges (see 
Fig.0. The out-of -plane spin structure factor meanwhile has 
enhanced correlations only at the spiral ordering wave vec- 
tors ±Q. These nontrivial properties distinguish the algebraic 
vortex liquid from other known spin liquids. Interestingly, the 
prominence of momenta Ki_2,3 in the theory appears to be 
consistent with recent series expansion studies of the Heisen- 
berg triangular antiferromagnet^^, which observe excitation 
energies at these wave vectors which are dramatically reduced 
relative to linear spin wave theory. Moreover, the prediction 
of "active" momenta ±Q and K1.2 in the anisotropic system 
seems to capture the neutron scattering data for Cs2CuCl4.i^ 

The phase diagram in the vicinity of the algebraic vor- 
tex liquid is also explored, and found to be rather rich. 
Nearby phases include the Kalmeyer-Laughlin chiral spin 
liquid, numerous magnetically ordered states including the 
coplanar spiral state, and variants of supersolids discussed 
recently3^2i^ 

In the range of anisotropy J' / J > 1.4, our treatment cap- 
tures the "square-lattice" Neel ordered state, which is the ex- 
pected ground state in this regime-^?i Here, we demonstrate in 
a particularly clear setting that anomalous "roton" minima in 
the excitation spectrum observed by series expansion studies 
can indeed be accounted for as low energy vortex-antivortex 
excitations.^''-''^ We further predict that these low-energy ro- 
tons may have still more dramatic effects in the easy-plane 
regime, which would be useful to explore using series expan- 
sions. 

The paper is organized as follows. The spin model and 
the dual vortex mapping are introduced in Sec. II. In Sec. 
Ill the fermionized vortex theory is developed. We first dis- 
cuss the "roton" excitations in the Neel phase arising when 
J'/ J > 1.4. We then obtain a low-energy effective theory for 
J'/ J < 1.4 which contains a description of the algebraic vor- 
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FIG. 1 : Triangular lattice and the dual honeycomb on which vortices 
reside. Spins shown illustrate a vortex. The spin exchange and vortex 
hopping amplitudes are generally anisotropic, with J' / J ~ t' /t. 



tex liquid. Sec. IV focuses on the properties of the algebraic 
vortex liquid, including its stability, symmetries, and dynamic 
spin correlations. The proximate phases of the algebraic vor- 
tex liquid are explored in Sec. V, and we conclude with a dis- 
cussion in Sec. VI. 



II. MODEL 

A. Easy-plane Spin Model 

We begin by considering an easy-plane, anisotropic spin- 
1/2 triangular antiferromagnet modeled by an XXZ Hamilto- 
nian with nearest-neighbor exchange, 

Ho = \Y. Jrr' [s+s;, + H.c] + J2 J^r'S^S^., (1) 

(rr'> (rr') 

where = S'Jf ± iS^ are the usual spin raising and lower- 
ing operators. As illustrated in Fig. ^ we take the in-plane 
exchange energy to be Jrr' ~ J along bold horizontal links 
and Jrr' = J' along diagonal links of the triangular lattice. 
The out-of -plane exchange is defined to be J^^, = 7 Jrr', with 
< 7 < 1 to satisfy the easy-plane condition. 

It is convenient to work with the easy-plane spin model 
recast in terms of quantum rotors by introducing an integer- 
valued boson number and its conjugate phase (^r- Upon 
identifying 5*^ — 1/2 and 5+ e"''"', the Hamiltonian 

reads 

Ho = ^ Jrr' COs{ipr - ipr') + U^{nr - 1/2)^ 
(rr') r 

+ ^ JrV(nr-l/2)(nr'-l/2). (2) 

(rr') 

The U term above energetically enforces the constraint of hav- 
ing either or 1 boson per site as required for modeUng a 
spin- 1/2 system. 
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The XXZ Hamiltonian Hq respects a number of internal 
and discrete lattice symmetries which we now enumerate. The 
model exhibits U( 1 ) spin symmetry and is invariant under time 
reversal T and a "particle-hole" transformation C which sends 
(S'^,S'y,S'^) ^ iS'^,-Sy,-S''). under these operations, 
the rotor fields transform as follows: 

U{1) : n^n, e"^ ^ e'('^+") (3) 
C : n^l-n, e*"^ ^ e"*"^ (4) 
T : n^l-n, e*"^ ^ -e"*^, (5) 

where a is a constant and T is an antiunitary operation which 
sends i — > —i. The model also preserves translations T^r by 
triangular lattice vectors 6r, as well as x-reflections TZx and 
inversions (i.e., tt rotations) Rt^ about a triangular lattice site. 
Rather than considering the x-reflections TZx, it will be useful 
for subsequent developments to work with a modified antiuni- 
tary reflection TZx = TZxCT. The latter operations transform 
the rotor fields as 

TZx-.Tir^nr', e*'^- ^-e*'^-', (6) 

where r' is an appropriately reflected coordinate. In the 
isotropic limit J = J', the XXZ Hamiltonian additionally pre- 
serves 7r/3 rotations i?7r/3 about a triangular lattice site; Rt^ is 
then no longer an independent symmetry since i^^yg = Rtt- 

In this paper we are interested in exploring the phase dia- 
gram accessible with the above XXZ spin model as a starting 
point. In particular, as the low spin and geometric frustra- 
tion strongly suppress the tendency to magnetically order, it 
is natural to ask whether spin-liquid phases can be realized 
with not too drastic perturbations to the model. To this end, 
we would like to derive an effective theory that governs the 
low-energy behavior of the spin system. This is, for instance, 
readily achieved for the integer-spin analogue of Eq. Q, and 
the phase diagram can be explored within a standard Landau 
analysis. However, for the spin- 1/2 system studied here a sim- 
ilar direct analysis of the spin model is hindered by the pres- 
ence of Berry phases. Consequently, obtaining a low-energy 
theory is largely intractable in this formulation. 

To proceed we utilize an alternative dual approach, in- 
troduced in the context of integer-spin systems in Ref. |25l 
wherein one considers a reformulation of the problem in terms 
of fermionized vortices. In this framework, the basic degrees 
of freedom one works with are vortices — topological defects 
in which the phases (fj- of the spins wind by 2tt around a trian- 
gular plaquette as in Fig. [2 — rather than the spins themselves. 
Although the vortices as defined are bosonic, it will prove ex- 
tremely useful to fermionize them in a manner famihar from 
the fractional quantum Hall effect via Chern-Simons flux at- 
tachment. Doing so enables us to obtain a low-energy dual 
theory, which as we will demonstrate leads naturally to a novel 
"critical" spin-liquid phase, the algebraic vortex liquid. 

B. Dual Vortex Mapping 

We proceed now to the dual vortex theory. We forgo the 
details of the duality mapping as these are provided in Sec. 



Ill of Ref. I25I in a very similar setting, and instead empha- 
size the important physical aspects of the dual theory. Im- 
plementing the duality transformation on the quantum rotor 
Hamiltonian Eq. (|2j,''' one obtains a theory of bosonic vor- 
tices with "electromagnetic" interactions hopping among sites 
of the dual honeycomb lattice depicted by the dashed lines in 
Fig-G] The vortices interact via a "vector potential" Oxx' G ^ 
and a conjugate "electric field" exx' which reside on honey- 
comb links and mediate a logarithmic vortex repulsion. Here 
X, x' denote nearest-neighbor honeycomb sites. (Throughout, 
we distinguish sites of the honeycomb and triangular lattices 
by the labels "x" and "r," respectively.) These dual gauge 
fields satisfy the commutation relation [exx'jOxx'] — i and 
commute on different links. In-plane spin components are en- 
coded in this formulation through the "electric field" and the 
vortices. The S'^ component of spin meanwhile appears as a 
dual "magnetic flux," 

S'; ^ -^(A X a)r, (7) 

ZTT 

where (A x a)r signifies a lattice curl of Axx' around the 
hexagon encircling site r of the triangular lattice. Although 
flxx' roams over the real numbers, the desired half-integer val- 
ues of 5*^ in Eq. Q are imposed energetically in the dual the- 
ory. 

In terms of a vortex number operator A'x and vortex cre- 
ation operator e'^", the dual vortex Hamiltonian can be ex- 
pressed as 

Hdual = - X! C0s(6'x - 6'x' " Oxx' " ftxx' ) (8) 

(xx'> 

together with a Gauss's law constraint for the "electric field", 

(A • e)x = iVx - 1/2. (9) 

Here a^^, is a static gauge field satisfying (A x aP)^ = tt, 
and (A • e)x denotes a lattice divergence of exx' at site x. 
Moreover, Ha describes the gauge field dynamics, 

na = ^xx'eL' +^^(A X a)l 

(xx') r 

+ ^ Jrr'(Axa)r(Axa)r', (10) 
(rr') 

withU ^ U/{2tt)\ J^%, = J^^,/{2tt)\ and J^^, = 27t^J' 
on the bold zigzag honeycomb links in Fig. ^ while j/xx' = 
2tt^J on vertical honeycomb links. 

The cosine term in Eq. (|SJl describes nearest-neighbor vor- 
tex hopping in an average background of tt flux per hexagon. 
This background "magnetic flux" is provided by the static 
gauge field a^x' and arises because is half-integer valued 
in the original spin model. (The average background flux for 
an integer spin system, in contrast, is trivial. ^^) The hopping 
amplitudes txx' are chosen to be anisotropic to reflect the spin 
exchange anisotropy. In particular, as illustrated in Fig.[2we 
take txx' = t on the bold zigzag links of the honeycomb and 
^xx' — t' on vertical links, with t' /t ~ J'/ J since vortices 
hop more easily across weak spin links than strong spin links. 
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TABLE I: Transformation properties of fields in the dual bosonic- 
vortex formulation under the discrete microscopic symmetries. The 
lattice coordinates, which also transform appropriately under the lat- 
tice symmetries, have been suppressed on all fields for brevity. 



Tst, R-k, Rtt/3 (isotropic limit) 


C 


7^.,T 


a" ^ a" 


a" ^ -a' 


a" ^ -a" 


a ^ a 


a —a 


a ~* —a 


e ^ e 


e — e 


e ^ e 


9^6 


-9 


9^-9 


N ^ N 




N ^ N 



An important feature of the dual theory is that with our 
conventions^^ the bosonic vortices are at half-filling, which is 
a direct consequence of the underlying frustration in the spin 
model. For example, in the classical VS x \/3 spin-ordered 
state we define the vortex number to be one on "up" triangles 
and zero on "down" triangles (or vice versa, depending on the 
chirality). The half-filling of the vortices becomes manifest 
upon reexpressing the dual Hamiltonian in terms of an uncon- 
strained electric field as follows, 

Tidual = Hdual + 5I(^x - l/2)14x'(A^x' ~ 1/2), (11) 
xx' 

where V^x' encodes the logarithmic vortex repulsion. Equa- 
tion M \\ clearly exhibits a vortex particle-hole symmetry. 

The transformation properties of the dual fields under the 
discrete microscopic symmetries can be straightforwardly de- 
duced as discussed in Ref. 25. These are summarized in Ta- 
ble U The continuous U(l) spin symmetry, which reflects 
conservation of 5*^, is not directly manifest in this formula- 
tion and is instead replaced by a conservation of dual gauge 
flux, (A X a). Additionally, the dual Hamiltonian has a U(l) 
gauge redundancy, being invariant under (oxx' + ^xx') ^ 
(flxx' + flxx' ) + — Ax' and (?x ^ ^^x + Ax for arbitrary 
Ax e M. 



III. FERMIONIZED- VORTEX FORMULATION 

Due to the finite vortex density together with the strong 
vortex interactions, the dual theory as it stands appears as in- 
tractable as the original spin model. There is, however, an 
important distinction between the dual vortex formulation and 
the original hard-core boson representation of the spin model 
that we can exploit. In the dual theory the vortices move in 
the presence of a dynamical gauge field which encodes the 
motion of the hard-core bosons. Thus, the dual vortex the- 
ory is in some sense a two-fluid model that describes simul- 
taneously both the dynamics of the hard-core bosons and the 
vortices. As such, it is possible to imagine a vortex moving 
together with a cloud of dual gauge flux, (A x a), which can 
in effect modify the statistics of the vortex-flux composite. In- 
deed, if the flux has strength ±27r, the composite particle will 
behave like a fermion due to the Aharonov-Bohm phase ac- 
quired from the dual flux under an exchange process. What we 
imagine is that the motion of the vortices and the dual flux can 



become dynamically correlated in such a fashion to be well 
represented (on intermediate length scales) by the dynamics 
of fermionic vortex-flux composites moving in the presence 
of the remaining dynamical gauge flux. 

This physical picture can be implemented by splitting the 
gauge flux into two pieces, a a + A, and attaching 2ti flux 
of (A X A) to the vortices with the help of a Chern-Simons 
term for A. An unfortunate but apparently unavoidable feature 
of this Chern-Simons approach is that we have to choose the 
sign of the attached flux, say, +2n rather than — 27r. One could 
contemplate an alternate formulation wherein the sign of the 
attached Chern-Simons flux is itself a dynamically fluctuating 
field, but we do not attempt to do so in this paper. 

Before proceeding to the details, we pause to comment on 
the usefulness and limitations of this approach. Working with 
fermionized vortices is expected to be legitimate for describ- 
ing physics in regimes where the vortex exchange statistics is 
unimportant. Consider, for instance, "insulating" phases of 
the vortices, examples of which include vortex crystals and 
"valence bond solids". Such phases were explored in integer- 
spin systems in Ref. 25 , and shown to coiTespond to magnet- 
ically ordered spin states. At the lowest energy scales, vortex 
density fluctuations are entirely frozen out, rendering their ex- 
change statistics unimportant; whether the vortices are treated 
as bosonic or fermionic is presumably inconsequential. 

On the other hand, once approximations are made to derive 
a low energy effective theory as we will do below, the vortex 
fermionization approach is expected to be least reliable when 
describing "vortex condensates". (Such vortex condensates 
correspond to paramagnetic states of the original spin sys- 
tem.) It is intuitively clear that trying to mimic the physics of 
Bose condensation will be challenging with fermionic fields, 
although this was the approach taken to describe anyon super- 
conductivity by a number of authors some years back. 

Here, we will be most interested in employing the fermion- 
ized vortex approach to access the "critical" algebraic vortex 
liquid. As we shall see, although the vortices are mobile in this 
phase, due to their long-range interactions vortex density fluc- 
tuations will be so strongly suppressed that the Chern-Simons 
flux attachment will be ineffective at modifying the behavior 
on long length scales. We will argue that the asymptotic low- 
energy physics of the algebraic vortex liquid is described by 
fermionic vortices minimally coupled to a gauge field medi- 
ating a long-range interaction (with Maxwell but no Chern- 
Simons term). This theory is usually referred to as QED3. 



A. Fermionization 

Formally, fermionization can be implemented by treating 
the vortices as hard-core bosons, replacing e**'= b]^ and 
A^x bl-bx ~ 0, 1, followed by a 2D Jordan- Wigner 
transformationjiS 

6], = 4exp[i ^ arg(x,x')A^x'], (12) 

iVx = 4^x = 4cfx. (13) 
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Here arg(x, x') denotes an angle formed by the vector x — x' 
with respect to an arbitrary fixed axis. 

The dual fermionized-vortex Hamiltonian takes the form 

(X1X2) 

+ R.c.]+na, (14) 
where we have introduced a Chern-Simons field 



A. 



x'7txi,X2 



[arg(x2,x')-arg(xi,x')]iVx' (15) 



which in Eq. il4\ resides on honeycomb links. Although we 
have included only nearest-neighbor hopping in the dual vor- 
tex Hamiltonian, one could also generic ally allow for small 
further-neighbor hopping terms allowed by symmetry. Upon 
fermionization, such terms similarly involve fermions coupled 
to a Chern-Simons field defined as in Eq. jl5> , but with xi and 
X2 further-neighbor sites. 

The transformation properties of the fermions and the 
Chern-Simons field can be deduced by examining Eqs. (I12> 
and (I15> . Tablellllsumm arizes the symmetry properties of all 
fields in this representation. According to Eq. (I12> . particle- 
hole symmetry sends d^.e^'''^, where for nearest- 
neighbor honeycomb sites Xi 2 the acquired phases satisfy 
7x1 - 7x2 = TT - 2(ylxiX2)- Here {A^,^^) denotes the mean- 
field value of the Chern-Simons field with (A'x) = 1/2 appro- 
priate for half-filled fermions. Since the Chern-Simons flux 
through a given hexagonal plaquette averages to 2tt, which is 
equivalent to zero flux, we take (A) = on nearest-neighbor 
links. Hence, in the table we implement particle-hole sym- 
metry by transforming d^ (— where j = 1,2 la- 
bels one of the two sublattices of the honeycomb. Note also 
that as discussed in Ref. |25l time reversal acts nonlocally on 
the fermions, and consequently we do not know how to faith- 
fully realize this symmetry in the continuum theory derived 
in the next subsection. In the last column of Table HIl we pro- 
vide a modified time reversal, Tfcrm, which acts locally on 
the fermion fields and corresponds to naive time reversal for 
fermions on a lattice. 



B. Vortex Mean Field and Low-energy Theory 

One advantage of working with fermionized vortices is that 
there is then a natural route to a low-energy effective the- 
ory. Namely, we start by considering a non-interacting "flux- 
smeared" mean-field state, ignoring fluctuations in the Chern- 
Simons and "electromagnetic" gauge fields and replacing the 
flux by an average background. Since the vortices are at 
half-filling, the Chern-Simons flux through each hexagon av- 
erages to 2tt, which is equivalent to zero flux on a lattice. 
Thus, the "flux-smeared" mean-field Hamiltonian describes 
free fermionic vortices hopping on the honeycomb in a back- 
ground of TT flux (due to flxx')" 




FIG. 2: Four-site unit cell chosen for the honeycomb. With our gauge 
choice, the static gauge field a^^^, is zero on the vertical links, while 
'^xx' ~ 7r/4 on the zigzag links directed along the arrows. The filled 
circle indicates our origin, which coincides with a triangular lattice 
site. 



Working out the fermionized-vortex band structure for 
Hmf exposes the important low-energy degrees of freedom in 
the mean-field theory. Doing so will enable us to derive a con- 
tinuum mean-field Hamiltonian, which will serve as the foun- 
dation on which we construct the full interacting low-energy 
theory by restoring fluctuations about the flux-smeared mean- 
field state. To this end, we diagonalize Hmf in momentum 
space assuming the four-site unit cell shown in Fig.|2l choos- 
ing a gauge with a'^^, — tt/A directed along the aiTows in the 
figure. Throughout we take as our origin a triangular lattice 
site denoted by the filled circle in Fig.|2l The band structure 
consists of four bands, two with positive energy and two with 
negative energy. Explicitly, the band energies at wave vector 
{kx, ky) in the Brillouin zone of Fig.|3lare given by 



2t^ ±2tJt'^ sin^ kx + + cos/c^ cos{V3ky)] 



Hmf — 



E 

(XX') 



H.c. 



(16) 



The spectrum has Dirac nodes at zero energy for t' /t < ^/2 
and is gapped for t' /t > \/2; we discuss these two cases sep- 
arately below. Our main focus will be on the former gap- 
less regime, where the algebraic vortex liquid arises. In the 
latter case with gapped vortices, which coiTesponds to the 
"square-lattice" Neel phase shown in Fig. |4] we will briefly 
discuss how within our flux-smeared mean-field treatment we 
can account for the anomalous "roton" minima in the excita- 
tion spectra observed by series expansion studiesi^iS 



1. t' /t > \/2.' Gapped vortices. 
"Rotons" in the frustrated square-lattice 

For t' /t > \/2, the half-filled fermionic vortices form a 
band insulator with the minimum band gap occurring at wave 
vectors Qi — (0,0) and Q2 = (tt, tt/a/S) in the Brillouin 
zone of Fig. |3l It is useful in this range of anisotropy to 
view the triangular system as a square lattice antiferromagnet 
with nearest-neighbor exchange J' and frustrating antiferro- 
magnetic exchange J along one diagonal direction as shown 
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TABLE II: Transformation properties of the fields in the dual fermionized-vortex representation. Symmetry properties of a ,a,e are the same 
as in TableU In the C column, j = 1,2 labels one of the two triangular sublattices of the honeycomb, and (A) refers to the mean-field value 
of Axx' with (A^x) ~ 1/2. The additional column Tfcrm corresponds to the naive time reversal for the lattice fermions and is not a symmetry 
of the vortex Hamiltonian. 





Tgr, Rtt, Rn/3 (isotropic limit) 




C 


T 






a" 


-a" 


-a" 


-a" 


-a" 


a 


a 


— a 


—a 


—a 


—a 


e — > 


e 


e 


— e 


e 


e 


dx 


d 


d 






d 


A^ 


A 


-A 


2{A) - A 


A 






FIG. 3: (a) Schematic energy dispersion with t' /t < \pl near the 
four Dirac points which compose the Fermi "surface" for the half- 
filled fermions in the flux-smeared mean-field state, (b) Locations 
of the Dirac points in the rectangular Brillouin zone corresponding 
to our unit cell choice in Fig.|2| The pairs of nodes Qi+, Qi- and 
Q2+, Q2- coalesce when i! jt — V2 and become gapped for t' /t > 
V2. 




AVL Neel order 
I \ ^ 

-1.4 t'/t~r/j 

FIG. 4: Phase diagram for the dual fermionized-vortex Hamilto- 
nian in the flux-smeared mean-field state. In the range of anisotropy 
t' /t > \/2, it is useful to view the system as a square lattice antiferro- 
magnet with nearest-neighbor exchange J' and frustrating coupling 
J along one of the diagonal links as shown above. 



in Fig. 13 The vortex insulator realized here corresponds to 
the square-lattice Neel state, which for sufficiently small J 
is the anticipated ground state. Since we are considering an 
easy-plane model, the spins order in the {S^, 5^) plane. The 
gapless Goldstone spin-wave at zero momentum is realized in 
the dual theory as a propagating "photon" mode in the elec- 
tromagnetic gauge fields. (The Goldstone mode at the order- 
ing wave vector present in a Heisenberg system acquires a 
gap in the easy-plane limit.) The Neel phase survives down 
to t' /t = V2, at which point the spectrum becomes gapless 
at Qi,2 signaling the destruction of the Neel order Series 
expansion studies^^ for the spatially anisotropic Heisenberg 
system find that the Neel phase survives in a similar range of 
anisotropy, J' /J > 1.4. 

Interestingly, excitation spectra for the Neel state calculated 
in series expansion studies of the Heisenberg system show 
significant deviations from spin- wave theory. ^^'^^ In terms of 
the standard square-lattice Brillouin zone notation, linear spin 
wave theory predicts identical excitation energies at momenta 
{tt/2, 7r/2) and (tt, 0) irrespective of the frustrating coupling 
J. Series expansions, on the other hand, obtain large en- 
ergy differences between these momenta due to a "roton" 
minimum in the spectrum at (tt, 0) which deepens as J in- 
creases (see Fig. 3 in Ref. 27). When J'/ J ~ 1.7, the ex- 
citation energy at (tt, 0) is roughly 27% lower than that at 
(^/2,7r/2)^ 

The anomalous minimum can be accounted for within our 
flux-smeared mean-field treatment as a low-energy vortex- 
antivortex excitation, thus substantiating the roton interpre- 
tation. Before proceeding, we want to note that Refs. Z/SJ 
consider — I excitations since the projection of the to- 
tal spin onto the Neel vector direction (assumed to be along 
S^) is conserved in the ordered phase of the Heisenberg sys- 
tem. There is no such spin quantum number in the easy-plane 
case, and we can characterize the excitations only by their mo- 
menta. If necessary, in the Heisenberg case there are always 
low-energy magnons near zero momentum and (tt, tt) which 
can be added to the vortex- anti vortex excitations discussed be- 
low to get the required spin quantum number 

To work more formally, consider the dynamical correlation 
of the operator at momentum q in the flux-smeared mean 
field theory. obtains contributions from vortex cuiTents 
whose circulation induces flux in the dual gauge field. The 
precise form of these vortex currents will be unimportant here, 
but can be obtained by constructing perturbations to the hop- 
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ping Hamiltonian which give rise to static gauge flux modu- 
lated at wave vector q. Genericafly, these contributions can be 
expressed as 

4 

E E>^&(k,q)4(kH(k-q). (17) 

a, 6=1 k 

Here a, b are band indices, k is summed over the Brillouin 
zone in Fig.lSjb), ^^(k) adds a fermion with momentum k 
in band a, and '^a,b are generally nonvanishing complex fac- 
tors. Consequently, the spin structure factor has contributions 
not only from spin-waves, but also from vortex-antivortex "ro- 
ton" excitations. In the ground state the two lower bands are 
filled, while the upper bands are empty. The excitation energy 
Arot (q) for a roton with momentum q is thus simply given 
by the minimum energy required to promote a fermion with 
arbitrary momentum k from an occupied band to a state with 
momentum k — q in an unoccupied band: 



1.5 



A,.ot(q) 



min{i;ompty 

k 



(k-q)-i;fmcd(k)} 



(18) 



This is straightforward to compute from the vortex band struc- 
ture. The result for Arot (q) along several cuts in momentum 
space is shown in Fig.|5]for three values of t' /t. To facilitate 
comparison with the series expansion results, we use square- 
lattice notation in the figure and display the same cuts as in 
Fig. 3 of Ref. 27, Note that unlike Ref. 27 which shows 
the lower edge for all excitations including spin waves, our 
figure shows only the vortex-antivortex excitations. We also 
point out that the vortex-antivortex excitation energies satisfy 
A,ot(q) = A,ot(q + (tt, 0)) = Ai.ot(q + (0, tt)), though this 
does not hold for general excitations. 

The most notable feature to observe in Fig.|5]is that the ro- 
ton excitation energy at (tt, 0) decreases as the frustration in- 
creases, consistent with series expansions, and becomes gap- 
less as the Neel order gets destroyed. The roton excitations 
at (0, 0) and (tt, tt) follow the same trend, though these low- 
energy rotons would be difficult to observe in a Heisenberg 
system due to the gapless spin-waves at these momenta. How- 
ever, since spin-waves at (tt, vr) are gapped in the easy-plane 
limit, significant deviations from spin-wave theory due to the 
roton at this wave vector are expected. Series expansions for 
an easy-plane system to search for this anomaly would be in- 
teresting, and could serve as a test for our explanation of the 
Heisenberg spectra. 



2. t' /t < \/2: Critical vortices 

For i = t' /t < V2 one finds that the Fermi "surface" 
for the half-filled fermionic vortices consists of four gapless, 
linearly dispersing Dirac points shown schematically in Fig. 
Ha). With our gauge choice these Dirac points occur at gen- 
erally incommensurate wave vectors Qi± and Q2± which can 
be written 



Qi± 

Q2± 



±{7r/2- 
±{7t/2- 



,0), 
,7r/V3) 



(19) 
(20) 
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FIG. 5: "Roton" excitation energy Arot(q) in the Neel phase along 
various cuts in momentum space. To facilitate comparison with Fig. 
3 from Ref. 21, we use the same square lattice wave vector nota- 
tion. As frustration increases, the roton excitation energy at (tt, 0) 
decreases, which is consistent with series expansion studies of the 



Heisenberg system,' 
order is destroyed. 



and eventually become gapless as the Neel 



with 



Q 



TT 

2 



(21) 



Fig-EIb) shows the positions of these wave vectors in the rect- 
angular Brillouin zone corresponding to our unit cell choice. 
The Dirac points have the familiar relativistic dispersion E w 



v^q^ + VyQ^, where q is measured relative to the nodal 

wave vectors Qu (L=l,2;l = +, — ). The velocities v^^y 
are in general anisotropic due to the anisotropy in the hopping 
amplitudes ^xx' and are given by 



1/2 



i4 



1/2 



(22) 
(23) 



Note some limiting cases: in the isotropic case i = 1 the ve- 
locities are equal; in the ID limit t ^ the spectrum becomes 
dispersionless in the y-direction; and finally, as i ^ \/2 and 
we approach the square lattice Neel state, the Dirac cones 
merge in pairs and flatten in the a;-direction, with Vx 0. 

For the purpose of exploring the low-energy physics of the 
theory, it suffices to focus only on low-energy excitations in 
the vicinity of the Dirac nodes. This can be achieved by ex- 
panding the fermion operators around the wave vectors Qu 
as follows. 



(24) 



Lla 



where x denotes sites of the honeycomb as before and £ = 
1, ... ,4 labels the corresponding site index in the unit cell 
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pictured in Fig.|2] On the right side of Eq. (I24> . are two- 
component spinors assumed to vary slowly on the lattice scale, 
with spinor components indexed by the label a Up 
to an overall uniform normalization factor, the "eigenvectors" 
^Lia ^^jj written 



$2-T 



e '12 



"-f 

s 


. 0. 

■ 0" 


s 

-t 

• 


-it 



s 
-it 


. . 



-if 
—s 


. 

-Q- 


s 
it 



(25) 



where s = 2 sin(Q/2). 

Using the expansion for the fermion operators in Eq. ( I24t . 
we obtain the following low-energy continuum description for 
the mean-field Hamiltonian, 



dx 'ijjli{-ivAcr'' - iVydya''^)ipLi, 



(26) 



where the flavor indices LI are implicitly summed and a^'^ 
are Pauli matrices that contract with the spinor indices [i.e., 
{a-'^4>)Lia = (^ap'^Liii]- Proceeding to the imaginai-y-time 
path integral formulation, the Euclidean Lagrangian density 
obtained from Eq. ( I26> can be written 



C-MF 



(27) 
(28) 



where the space-time index /i = 0,1,2 is defined so that 
do,i,2 = dr.x.y and we have rescaled the spatial coordinates 
to absorb the anisotropic velocities ti 



7' 



are given by 7° = 



7 



V ^,2 



a", 7 



The Dirac matrices 
(T^ and satisfy 



the usual algebra {7^,7"} = 2(5^"^. We will also frequently 
use Pauli matrices /i*^ and r*^ which contract with the flavor 
indices L and L i.e.. 



Ll 



(29) 



Upon resurrecting the vortex interactions and gauge field 
fluctuations about the mean-field state, we obtain the desired 
low-energy theory. 



^ = ^L/7^(9,,-«a 
1 

+ Ye 



iAMi 



i 

47r^ 



2 {<^fiu\d„axf + -^f-^j^uxA^d^Ax + Af ■ (30) 



Equation (13 0> describes four flavors of two-component Dirac 
fermions ipu (corresponding to the four Dirac points) min- 
imally coupled to a non-compact U(l) gauge field and a 
Chern-Simons field A^. The gauge field mediates the loga- 
rithmic vortex repulsion, while the Chern-Simons terms above 
enforce the flux attachment to the fermions, thereby restoring 
the original bosonic vortex exchange statistics. The form of 
the Maxwell term above is only schematic and ignores the un- 
derlying lattice anisotropy. Finally, Ca represents symmetry- 
allowed four-fermion terms arising from short-range parts of 
the vortex interactions in the microscopic model. We furnish 
an explicit form of in Sec.lVl 

Table |lll] displays the transformation properties of the 
continuum fermion fields under the microscopic symmetries 
(with Tform rather than T due to subtleties mentioned above). 
With out gauge choice for a^^' in Fig- 13 the two translations 
given in Table IIIII are realized as follows. The first, Ti, corre- 
sponds to a simple translation of fields by 5r = x. The sec- 
ond, T2, corresponds to translation by = — l/2x + V3/2y 
and must be accompanied by a gauge transformation. Specif- 
ically, T2 is realized by first implementing the required gauge 
transformation by sending 



(31) 
(32) 



where rix^y are integers labeling the unit cell to which site x 
belongs, and then translating as follows. 



d. 



x,l/2 
^x,3/4 



<:+i5r,2/l 
<:+(5r,4/3- 



(33) 
(34) 



Particle-hole symmetry and fermionic time reversal similarly 
require gauge transformations, as do rotations in the isotropic 
limit J = J'. 

Translations, reflections, and particle-hole symmetry pro- 
hibit all possible fermion bilinears from appearing in Eq. 
( I30> except iififi^T^a-y^ — n^r^a^^ and ifj^p = a^-tp. 
The first of these bilinears is a perturbation to the Hamilto- 
nian which modifies the ratio of vortex hopping amplitudes 
t' /t. This has the trivial effect of shifting the x-components 
of the nodal wave vectors Ql;, which are not protected in 
an anisotropic system. In the isotropic limit, itp t'^ tp is 
eliminated by rotation symmetry. Deducing the fate of ipip, 
which respects all symmetries except Tfcrm, requires more 
care and will be discussed in Sec. IIV Al We will argue 
that adding this term to the action drives the system into the 
Kalmeyer-Laughlin chiral spin-liquid, which breaks physical 
time-reversal symmetry. Thus should be excluded if we 
are to describe a time-reversal invariant state. 

At this point it is worth emphasizing that upon consider- 
ing the simplest flux-smeared state, we have already arrived 
at a mean-field description of the "critical" algebraic vortex 
liquid which is the main focus of this paper. The "critical" na- 
ture of this state follows from the gaplessness of the fermionic 
vortices, which in turn allows for gapless spin excitations as 
we will discuss below. Many properties of the AVL, such as 
the momentum-space locations of the low-energy spin excita- 
tions, can in fact be deduced from the mean-field theory. By 
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TABLE III: Transformation properties of the continuum fermion fields ijj. Symmetries T\ and T2 correspond to translations by (5r = x and 
5v = — l/2x + \/3/2y, respectively. Moreover, Tfcrm corresponds to the naive time reversal for fermions on the honeycomb rather than the 
physical spin time reversal. 





Ti 








c 


'^crm 


Rt^/'j, (isotropic limit) 



















studying the effective Lagrangian Eq. j30L we will attempt to 
go beyond such a mean-field analysis. In particular, in the fol- 
lowing section we will address the stability of the AVL when 
fluctuations about the flux-smeared state are incorporated and 
make quantitative predictions for various spin correlations in 
the AVL. Moreover, with this effective theory in hand we can 
also explore the phase diagram in the vicinity of the AVL. The 
algebraic vortex liquid has a number of interesting proximate 
phases, some of which we explore in Sec.fVl 

To summarize, the mean-field phase diagram along the spa- 
tial anisotropy axis t' /t is shown in Fig.0] We will focus on 
t' /t < \f2 for the remainder of the paper. 



IV. ALGEBRAIC VORTEX LIQUID 

A. Effective theory of the AVL— QED3 

We begin our detailed characterization of the algebraic vor- 
tex liquid by examining the continuum theory describing fluc- 
tuations about the critical flux-smeared mean-field state. The 
full interacting theory is described by the effective Lagrangian 
Eq. ( I30> . To ascertain response properties of the spin system, 
we add an external probing field A^^'^ which couples to the 
three-current of the hard-core bosons. In the dual vortex for- 
mulation prior to fermionization, this three-current is given by 
6j = (V X a)/2Ti. Upon fermionization, we introduced an ad- 
ditional Chern-Simons field, and attached flux (V x A) = 27r 
to the vortices. As discussed at the beginning of Section III, it 
is convenient to view the Chern-Simons field as being "part" 
of the original gauge field, that is a — > a = a + A. The idea 
is that it is the physical hard-core boson current which is be- 
coming correlated with the motion of the vortices, forming a 
vortex-flux composite which behaves as a fermion. Based on 
this picture, it is reasonable to assume that the Chern-Simons 
gauge flux carries the U{\) charge of the hard-core bosons, 
and to couple in the external field via 

/:cxt = • [(V X a) + (1 - k)(V X A)], (35) 

with K = 0. Traditional application of Chern-Simons the- 
ory would take instead k = 1, which corresponds to the as- 
sumption that the Chern-Simons flux is "fictitious" rather than 
physical and hence carries no quantum numbers. As we shall 
see below, the choice k = is preferable, being essentially 
equivalent to replacing physical time reversal invariance by 
Tfcrm- But for now we keep k, as an arbitrary parameter. 



Since the Dirac fermions in the low-energy continuum the- 
ory given in Eq. ( I30> couple only to = + A^, it is in- 
structive to rewrite the Lagrangian in terms of this sum field. 
The full Lagrangian can then be cast in the following appeal- 
ing form: 

^ — -CqEDS + ^CS + -Cint + ^Gxt, (36) 

with 

+A/ , (37) 

and 

£cs = ^A- (V X A). (38) 

47r 

Here, £qed3 describes non-compact quantum electrodynam- 
ics in 2 + 1 dimensions (QED3) with = 4 flavors, which is 
coupled to the Chern-Simons Lagrangian by an interaction, 

Ant = -4(Vxa).(Vx^). (39) 
The external probing field takes the form 

^cxt = -^A'^^' • (V X 5) + n^A-'^' ■ (V x A). (40) 

ZTT ZTT 

Notice that with the choice k = 0, the external source field 
only couples to d. 

Before discussing the effects of the interaction term, we 
briefly review the behavior of QED3, which has been widely 
studied in a variety of contexts.^^i^^i^'^i^^-^''-^^-^'^-^'^-^" The fixed 
point with = in which gauge fluctuations are entirely 
suppressed and the fermions are essentially free is unstable, 
so that QED3 is inherently a strongly interacting field theory. 
Consequently, to make progress analytically one must modify 
the theory in a manner which provides a controlled limit. An 
often used approach is the large-iV limit, where one general- 
izes to a large number of fermion flavors. Starting from the 
infinite- limit, one can then perform a controlled analysis by 
perturbing in powers of 1 /N. This approach can be cast in the 
form of a renormalization group treatment, and an important 
feature is that the gauge field a scales like an inverse length 
and due to gauge invariance does not pick up an anomalous 
dimension. The scaling dimensions of the symmetry-allowed 
four-fermion interactions in Eq. (|^} do generally acquire an 
anomalous dimension, which can be computed perturbatively 
in inverse powers of N. For large enough N all four-fermion 
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terms are irrelevant, and QED3 thus realizes a nontrivial sta- 
ble critical phase. For N < N^., with some unknown Nc, it 
is believed that four-fermion terms become relevant, leading 
to spontaneous fermion mass generation and the destruction of 
criticality (except with fine-tuning). While Nc can in principle 
be deduced in QED3 simulations, recent studies are inconclu- 
sive as to whether Nc lies above or below the = 4 case of 
interest in the present work.^'' '*" We will assume henceforth 
that Nc < 4. Further numerical simulations and higher-order 
calculations in would be useful for justifying (or negat- 
ing) this assumption. 

We are now in position to consider the effects of the cou- 
pling dnt between QED3 and the Chern-Simons Lagrangian. 
Since Ccs is Gaussian, it can be viewed as the fixed point of a 
renormalization group transformation in which A is rescaled 
like an inverse length. For large e, the effects of the interac- 
tion can then be studied perturbatively. Since £int is quadratic 
in the gauge fields and involves two derivatives it has scal- 
ing dimension 4, and is formally irrelevant. The fixed points 
described by QED3 and Chern-Simons theory evidently de- 
couple at low energies. The physics here is that due to the 
logarithmic vortex repulsion which strongly suppresses vor- 
tex density fluctuations, exchange statistics play only a minor 
role at criticality. In the next subsection we will employ QED3 
to access the properties of the critical algebraic vortex liquid 
phase. 

But caution is necessary. For £qed3 fermionic time rever- 
sal symmetry precludes the generation of a fermionic mass 
term ipip, which respects all symmetries of the Lagrangian 
except Tfcrm- The Chern-Simons Lagrangian Ccs^ however, 
is not invariant under Ifoim- As a result, once the two the- 
ories are coupled, despite the irrelevance of this coupling a 
small fermion mass term will presumably be generated, be- 
ing no longer symmetry-protected. Tracing back its origin, 
we see that the sign of the vortex mass will be determined by 
the sign of the flux that was attached upon vortex fermioniza- 
tion. As we will discuss in Sec. IV A| this mass term drives the 
system into the Kalmeyer-Laughlin chiral spin-liquid, which 
breaks physical spin time-reversal symmetry. Thus, in or- 
der to correctly implement a renormalization group analysis 
that faithfully respects the physical time reversal symmetry 
of the original spin model, we must maintain masslessness of 
the fermions. We will proceed under the assumption that the 
physically correct procedure is to tune a small bare mass term 
to cancel the effects of the irrelevant coupling as it scales to 
zero — that is, to tune the fully renormalized mass term to zero. 
The resulting massless and critical QED3 gives us a descrip- 
tion of the time reversal invariant algebraic vortex liquid. 

Subtleties associated with time reversal invariance are also 
apparent in the Hall conductivity of the original hard-core 
bosons, which we now briefly discuss. Once the Chern- 
Simons Lagrangian Cqs has decoupled, one can readily per- 
form the Gaussian integration over A, which gives, 

= — ^-A^'^t . (V X 5) - K^—A"""^ ■ (V X A"""^). (41) 
27r 47r 

This form shows that the conductivity tensor of the original 



hard-core bosons, aap, is given by 

fv 9 

^"^=(2^+2^'"^' ^^^^ 

where p^^ is the resistivity tensor for the fermionic vortices 
described by QED3 and e is the antisymmetric tensor. Since 
the Hall resistivity for the fermions vanishes in QED3 due to 
^crm symmetry, the Hall conductivity for the bosons in the 
critical AVL phase is given by (j^y — K'^/27r. Recall that the 
parameter k gives a measure of how much bosonic charge is 
ascribed to the statistical flux attached to the vortices. 

Generally, the Hall conductivity is not a low-energy prop- 
erty of a physical system, and as such can be non-universal 
even at a critical point or in a critical phase. However, since 
the original spin model is time reversal invariant, the Hall con- 
ductivity must vanish, at least in the absence of any sponta- 
neous symmetry breaking. Notice that the required vanishing 
of the Hall conductivity follows provided we take the param- 
eter K — 0. As discussed above and in Section 111, the choice 
K = corresponds to assuming that the Chern-Simons flux 
attached to the vortices carries a non-vanishing boson charge. 
On the other hand, if the statistical flux is presumed to carry no 
charge, one has k = 1 and a non-vanishing Hall conductivity, 
with a sign set by the sign of the attached statistical flux. The 
former choice, k — 0, gives us a way to access a time reversal 
invariant state with zero Hall conductivity independent of the 
sign of attached flux. Physically, with the fermions coupling 
to the sum a = a + A, the motion of the vortices leads to a 
"screening" of the attached flux V x ^, by the fluctuating flux 
V X a. At long wavelengths the total flux surrounding each 
vortex, V X a, which is proportional to the full boson current 
when K = 0, vanishes. The vortices are charge neutral and the 
Hall effect vanishes. 

It would clearly be desirable to have a method for fermion- 
ization involving flux attachment in a more democratic fash- 
ion which treats +27t and —2tt in an exactly equivalent man- 
ner. But in the absence of such an approach, we must content 
ourselves with using £qed3 together with the assumption of 
K = to describe the properties of the time-reversal invariant 
AVL phase. 

B. Symmetries of the AVL 

The critical QED3 theory proposed to describe the AVL re- 
spects all symmetries in Table lllll and also has a dual global 
U(l) symmetry under tp e'"-!/" reflecting conservation of 
vorticity. Due to the assumed irrelevance of four-fermion 
terms in the scaling limit, the theory also possesses an emer- 
gent global SU(4) flavor symmetry, being invariant under ar- 
bitrary SU(4) flavor rotations of the form Ut/j, with 
— U^^. The 16 conserved three-currents associated with 
the U(l) and SU(4) symmetries can be compactly written 

^a;3 = V^7>"r'3^ (43) 

and satisfy d^J'^p = 0. Here the indices a,/? range from 
to 3; /x" and r° are identity matrices; and ji^ and are 
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FIG. 6: Momenta carried by the fermion bilinears in Table llVI whose 
correlations are enhanced by gauge fluctuations at the AVL fixed 
point. 



Pauli matrices defined as in Eq. ( I29> . The U(l) conserved 
current is Jqq, while the remaining 15 currents constitute the 
SU(4) conserved currents. The 48 fermion bilinears compris- 
ing are prohibited from acquiring an anomalous dimen- 
sion. We will be primarily interested in the remaining 16 
bilinears, whose correlations are enhanced by gauge fluctu- 
ations; their transformation properties are supplied in Table 
II VI Figure |6l displays the set of momenta caiTied by these 
enhanced fermion bilinears, which coiTespond to the leading 
gapless vortex-antivortex excitations. The wave vectors in the 
figure are explicitly given by ±Q, where 

Q = (2(5 + TT, 0) (44) 

and Q is defined in Eq. i21\ . Ki 2 = (tt, =F7r/v^), K3 = 
(0,27r/V3), and ±Pj = ±(Q + K^). We will often refer 
to ±Q as spiral ordering wave vectors, since in the isotropic 
limit these correspond to the \/3 x \/3 order. 



C. Dynamical spin correlations in the AVL 

We turn now to the dynamical spin correlations in the AVL. 
Due to the gaplessness of the fermionic vortices, the AVL ad- 
mits universal power-law coiTelations in the spin structure fac- 
tor. To extract these spin correlations from our dual theory, we 
need to first identify the operators in QED3 which correspond 
to and 5+. We discuss the correlations of and 5+ sep- 
arately below. 



carry the same quantum numbers. More precisely, due to sub- 
tleties with realizing physical time-reversal in QED3, we re- 
quire that contributing fermion bilinears have the same quan- 
tum numbers as V x a. These bilinears arise microscopically 
from gapless vortex currents that induce gauge flux modula- 
tions at finite wave vector, which are analogous to the "rotons" 
discussed in Sec. llllBTI 

From such an analysis, the continuum expression for 
takes the form 



5" 



A X a 
2n 



ss ■ 



H.c] 



(45) 



where Aiss and Aigg are enhanced fermion bilinears from 
Table HVl that carry momenta ±Q. It can be readily verified 
using the table that the right-hand-side has the desired symme- 
try properties. The ellipsis in Eq. (145 > represents terms aris- 
ing from non-enhanced fermion bilinears (i.e. , those which are 
part of the SU(4) conserved currents) and higher-order contri- 
butions. For completeness, we note that these non-enhanced 
bilinears carry momenta Ki 2,3 and ±Pi_2.3 in Fig.|6l Thus 
the coiTelations at momenta Kj and ±Pj have scaling di- 



mension An 



2. 



Due to enhancement from gauge field fluctuations, the 
fermion bilinear A^gg in fact provides the dominant power 
law in the correlations. Near momenta ±Q the spin 
structure factor scales as 



5^^(k = ±Q + q,t^) 



(^2 _ q2)l-r,„„„/2- 



(46) 



The anomalous dimension rjanh is that of an enhanced fermion 
bilinear in QED3 and can be estimated from the leading 1/iV 
result-21 



Vcnh 



128 



(47) 



Setting = 4 yields ?7cnh ~ 1-92, and a scaling dimension 
Acnh = (1 + '?cnh)/2 ~ 1.46. At all other wave vectors 5""^ 
exhibits subdominant power laws. 

The fact that the leading correlations occur at momenta 
±Q suggests proximity of "supersolids" to the AVL. Such 
states are characterized by concurrent S"^ and 5+ order, and 
are nontrivial in an easy-plane system where typically only 
in-plane spin order occurs. As we will show in Sec. IV B 21 
supersolid phases indeed emerge naturally out of the AVL. 



S correlators 



1. correlators 

From our microscopic identification in Eq. 0, it is clear 
that near zero momentum appears in QED3 as the con- 
served dual gauge flux, (V x a) /2tt. Since a is massless in the 
critical theory, exhibits power-law correlations at zero mo- 
mentum with scaling dimension 2. At other wave vectors, 
receives contributions from fermion bilinears in QED3 which 



Since the spin raising and lowering operators S add spin 
±1, the corresponding operators in QED3 are "monopole in- 
sertions" which add ±27r gauge flux. Our goal here will be 
to construct a continuum expression for 5+ in terms of these 
monopoles and from this extract the leading in-plane spin cor- 
relations, much as we did for above. Monopole operators 
in QED3 were discussed in a very similar setting in Ref. 13 
and here we shall only highlight the main points. We will as- 
sume that the added ±27r flux is spread smoothly over a large 
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TABLE IV: Transformation properties of tiie 16 bilinears whose correlations are enhanced by gauge fluctuations at the AVL fixed point (we do 
not show separately A4^gg and Vl 2 3)- In the second column, Q is the spiral ordering wave vector defined in Eq. <44> . Kj lie on the midpoints 
of the Brillouin zone edges, and Pj = Q + Kj . Figure|6|displays the set of momenta carried by the enhanced fermion bilinears. 
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area compared to the lattice unit cell, and treat the flux as a 
static background. This flux alters the fermionic spectrum, 
and in particular gives rise to four zero-energy modes, one for 
each fermion flavor in the continuum. The zero-mode wave 
functions can be obtained by first modifying the mean-field 
Hamiltonian density in Eq. (I26> as follows. 



HMF.g ~ -i^Uid^ - iaDa^ + {dy - ial)<jy]^Li- (48) 



Here a| y is the vector potential giving rise to 2'n:q flux, with 
q — ±.1 the monopole "charge". Focusing only on the zero- 
modes, we then replace 



modes: 



(49) 



where (pu.q are the desired zero-mode wave functions and 
the operator fLi,q annihilates the coiTesponding zero-mode. 
Choosing the Coulomb gauge for a% y and assuming an az- 
imuthally symmetric flux distribution centered around the ori- 
gin, it is straightforward to show that the zero-mode wave 
functions araii 



|x| (0, 
Ixlll 



(50) 
(51) 



The transformation properties of the zero-mode operators 
fLi,q can be deduced from Eq. ( I49> and the transformation 
properties of V'l/ in Table HiH the results are shown in Ta- 
ble|Vl 

Since the fermions are at half-filling, physical {i.e., gauge- 
invariant) states must have two of the four zero-modes oc- 
cupied. Thus we need to consider six distinct monopole in- 
sertions. It will be useful to define the following translation 
eigenoperators which add fermions to two of the four zero- 
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= fl 
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^L.q 


= fl 









(52) 



For convenience, the transformation properties of these oper- 
ators are also displayed in Table IVl 

We now introduce the monopole operators by specifying 
their action on the ground state with no added flux, denoted 
|0). First, we define monopole creation operators A/j, which 
insert +27r flux and fill the zero-modes as follows. 



AfjiO) 



e"^^Fl^^\DS,+l) 



R/L 



|0) = e»«/-F^/^^^Ji?5,+l 



(53) 
(54) 
(55) 



Here j runs from 1 to 3 and \DS, +1) is the filled negative- 
energy Dirac sea with +27r flux inserted and all zero modes 
vacant. The Hermitian conjugate operators are required 
to add the opposite momentum and flux to the ground state. 



AfolO) 
MM 



= e^P^Fl_,\DS,-l) 
= e"'^Fl_^\DS,-\) 



L/fl,-l 



\DS,-l) 



(56) 
(57) 
(58) 



where \DS^ —) is the filled negative-energy Dirac sea with 
—2tt flux inserted. It is important to note that the phases a 
and [3 in the definitions above are arbitrary, and can be spec- 
ified to our convenience so as to construct operators with the 
desired transformation properties. 
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TABLE V: Transformation properties of tlie zero-mode operators fu.q in tlie ciiarge q = ±1 monopole sectors. Tiie quoted transformations 
were obtained by employing the Coulomb gauge for the added ±2n gauge flux. Also shown are the transformation properties of the six 
translation eigenoperators , defined in Eq. <52> which add fermions to two of the four zero-modes. 
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To obtain a continuum expression for in terms of the 
monopole operators, we need to determine the quantum num- 
bers they carry. This in turn requires knowing the transforma- 
tion properties of \DS, (?) (0|. We take up this rather involved 
issue in Appendix|X| By employing general relations among 
the symmetries (such as = 1, etc.), we first establish that 

ri,2 : \DS,q){0\^~\DS,q){Q\ 

: \DS,q){0\^iUDS,q){0\ 

: \DS,q){0\^C^\DS,q){0\ 

C : \DS,q){0\^Y[f\^_^\DS,~q){0\ 

Aa 

Tferm : \DS,q){0\^\DS,-q){0\ 

R^/s : \DS,q){0\^C^e^'i^^''\DS,q){0\. (59) 

Here and are signs which in principle are fixed, but can 
not be determined using only symmetry relations. The speci- 
fied action under rotation applies only in the isotropic limit. 

This information is sufficient to determine the momenta 
carried by the monopoles: AIq caiTies zero momentum, AIj 
carries momentum Kj on the midpoint of a Brillouin zone 
edge (cf. Fig. Q, and Mh/l caiTy momenta at the spiral or- 
dering wave vectors ±Q. For other symmetries, however, a 
more careful analysis is required. In fact, by examining the 
symmetry of the monopoles under inversion, one can show 
that it is impossible for all six bare monopole operators to 
contribute to 5*+. Under inversion, we have A/q Cn^Ia 
while Mj —C,TrMj. But the Fourier components of 5+ 
at zero momentum and Kj are all even under inversion. De- 
pending on the sign (t^, either the bare monopole Mq or the 
three monopoles Alj must therefore be excluded on symmetry 
grounds from a continuum expression for 5*+. (Since Mn/j^ 
are not diagonal under inversion, one can always choose the 
phases aji/L and (3jif]^ to construct operators that transform 
like at momenta ±Q). 

To determine the remaining ambiguities we appeal to nu- 
merical studies of monopole insertions. Specifically, we diag- 
onalize the mean-field hopping Hamiltonian on a finite system 
with arbitrary flux insertions to obtain the single-particle ener- 
gies and wave functions. With these wave functions in hand, 
it is then possible to obtain the inversion and reflection prop- 
erties of \DS,q){Q\. These numerics are discussed in more 



detail in Appendix |X| To ensure geometry-independence of 
the results, a variety of boundary conditions and system sizes 
were considered. In all cases, we find that (^^r = = ^ 1- In 
particular, we conclude that Mq by itself does not contribute 
to a continuum expression for S^. 

We now have enough information to determine unambigu- 
ously all quantum numbers caiTied by the monopoles. It is 
straightforward to show that the phases appearing in Eqs. i53i 
through j58> can be chosen so that the monopoles transform 
as shown in Table lVll (Note that under physical time-reversal 
5"+ — ^ ~S^, whereas Tfcrm sends +Ma. This is 

not too surprising, however, given that the U(l) spin symme- 
try is not manifest in the dual theory.) The desired continuum 
expression for 5+ can then be written as follows. 



5-+ ^ [e-'^-'^Alj^ + e'^-^All] 



,(60) 



where the ellipsis represents subdominant contributions. The 
momenta ±Q and Kj carried by the monopoles on the right- 
hand-side are sketched in Fig.0 

Monopole operators are known to have nontrivial power- 
law correlations in large-iV QED3, each with identical scal- 
ing dimension A,n ~ 0.26iV.'*^ This fact leads us to a non- 
trivial prediction for the in-plane spin structure factor 5^ in 
the AVL. Namely, 5*^ exhibits the same universal power-law 
correlations at each of the five momenta Uj of Fig. Q This re- 
markable property stems from the enlarged global SU(4) fla- 
vor symmetry enjoyed by the AVL. For wave vectors near Ilj, 
the in-plane structure factor thus scales as 

^+-(k = n, + q,.) ^ (61) 

The anomalous dimension t]^ is given in the large-A^ limit by 



Vn 



0.53A^- 1. 



(62) 



Setting iV = 4 yields r^m w 1.12 and A„i w 1.04. While the 
scaling dimension is the same at each wave vector, it is im- 
portant to keep in mind that the amplitudes Aj can vary sig- 
nificantly at different momenta. For example, near the limit 
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FIG. 7: Momenta at which the dynamical spin stracture factor S'^ 
exhibits dominant power law correlations with the same exponent 
in the AVL. The leading power law in the structure factor, by 
contrast, occurs only at momenta ±Q. 



of decoupled chains the amplitude at K3 should be much sup- 
pressed relative to the other four wave vectors, which are near 
kx = TT where most of the activity would be expected. 

We note here that while exclusion of the bare monopole Mq 
from a continuum expression for 5*+ was not obvious at the 
outset, this conclusion is quite reasonable physically in light 
of the spin coiTelations discussed above. If this exclusion did 
not occur, then the dynamic spin structure factor 5*^ would 
exhibit the same power-law coiTelations at zero momentum 
and the five wave vectors in Fig.Q This would be quite sur- 
prising given that one would intuitively expect subdominant 
correlations at zero momentum in an antiferromagnet. 

The locations of the leading in-plane coiTelations in the 
AVL are suggestive of proximity to magnetically ordered 
phases involving condensation of 5+ at the momenta shown 
in Fig.0 We will explore some of these states below. 



actions are nevertheless expected to generate fermion masses 
and destroy the AVL. We will focus, in particular, on states 
arising from the generation of enhanced fermion mass terms, 
as the AVL will likely be more susceptible to realizing such 
states. Moreover, for simplicity we will consider the spatially 
isotropic limit J = J'. In this case, the four-fermion terms 
£4f can be written in terms of enhanced bilinears in Table II Vl 
as follows, 

Ui = uiM]iL + u2M^^^^ + u:iM\sMss 
3 

+ Y^^UiK^+u^Kf + U(iV]V,) + Cii,r... (63) 

The last term £4f ,nc represents four-fermion interactions com- 
posed entirely of non-enhanced bilinears, which will not be of 
interest here. Though not unique, this form provides a useful 
organization of the four-fermion interactions based on their 
translation and rotation properties. Our exploration below is 
by no means intended to be exhaustive; rather, our aim is to il- 
lustrate some representative examples of the proximate phases 
that can be analyzed in the fermionized vortex theory. 



A. Kalmeyer-Laughlin spin liquid 

Consider first the addition of a mass term niAiKL = rmpip, 
which is favored by a large negative ui interaction above. This 
mass respects all symmetries except Tfcrm, and drives the sys- 
tem into a v — 1/2 fractional quantum Hall state for the origi- 
nal bosons, which breaks physical time-reversal symmetry. In 
other words, the physical spin state obtained by the addition 
of mAiKL is the Kalmeyer-Laughlin chiral spin-liquid.^ *^ 

To demonstrate this, let's first integrate out the massive 
fermions. Since all flavors have the same mass m with the 
same sign, integrating out the fermions induces a Chern- 
Simons term for (a + A)^. The Lagrangian is then 



V. PROXIMATE PHASES TO THE ALGEBRAIC VORTEX 
LIQUID 

In this section we explore the neighboring phases of the 
AVL that are encoded by the effective Lagrangian ( I30> . Re- 
sponse properties of the bosonic spin system will be obtained 
by introducing an external probing field A°^^ that couples in 
the dual theory via Eq. ( I35> . As discussed in Sec. llV Al we 
will assume /t = so that A°^^ couples to both the original 
boson currents Sj = (V x a)/27r and the Chern-Simons flux. 
This exploration will provide some guidance as to where in 
the phase diagram the AVL lies, and is also useful because a 
study of the phase diagram within a direct Landau analysis of 
the spin model is hindered by Berry phases. 

Descendants of the AVL are obtained by giving the 
fermions a mass, which destabilizes the vortex liquid lead- 
ing to a wealth of possible states. Here, we will restrict our 
attention to nearby states favored by the interactions . Al- 
though we postulated above that such terms are irrelevant in 
the AVL critical theory, sufficiently strong four-fermion inter- 



C-a,A 



_L(Vxa)^ + ^A.(VxA) 



isign(m) 



27r 



27r 



(a + ^) • [V X (a + A)] 
[V X (a + A)] . 



(64) 



The spectrum for the above Lagrangian is gapped, which can 
be verified by integrating out the Chern-Simons field A. Inte- 
grating out further the gauge field a, we arrive at an effective 
Lagrangian for the probing field 



£eff = -sign(m)^A-* • (V x A^*), 



(65) 



with (Jxy = vj^-K = I/Att. Thus Eq. i65i characterizes the 
response for a i/ = 1/2 fractional quantum Hall state of the 
original bosons as claimed. We remark here that this physics 
results with either sign for the mass ni. Had we instead chosen 
to couple ^4°''* only to the original boson currents [i.e., with 
K = 1 in Eq. ( I35H . then a particular sign of the mass would 



15 



TABLE VI: Transformation properties of the six monopole operators in QED3. 
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have to be chosen relative to the sign of the Chern-Simons 
flux in order to recover the Kalmeyer-Laughlin state. Once 
again, we see that endowing the Chern-Simons flux with bo- 
son charge leads to response properties of the spin system that 
are insensitive to the direction of flux attachment as desired. 

What is the nature of the gapped excitations in this phase? 
Consider acting on the ground state with the fermion field . 
The added fermion couples to the sum a = a + A. By examin- 
ing the action obtained by retaining d and integrating out the 
Chern-Simons field, we see that the system dynamics binds 
Ax a — — TT flux to the fermion (which also carries 2tt Chern- 
Simons flux). Thus, the fermion is turned into a semionic exci- 
tation carrying spin- 1/2. This is precisely the gapped semionic 
spinon in the Kalmeyer-Laughlin state. 



B. Magnetically ordered phases 

The remaining states we consider arise from generating 
specific fermion mass terms of the form mipWip, where W 
has two +1 and two —1 eigenvalues. Hence, two fermion 
modes have mass m, while the other two have mass — m. In 
all such phases, the vortices are "insulating," and the "pho- 
ton" in the dual gauge field a can freely propagate. The 
gapless photon is revealed upon integrating out the massive 
fermions, which induces only a generic Maxwell term for 
the field d — a + A. These vortex insulators correspond to 
magnetically ordered phases of the spin system. The gapless 
photon is the Goldstone spin-wave at zero-momentum arising 
from the broken continuous U(l) spin symmetry. Moreover, 
the probing field A''^^ is massive here, which is the "Meiss- 
ner effect" expected for the superfluid phase of the original 
bosons. Our objective below will be to disentangle the spin 
order that arises in different vortex insulators. As we will see, 
magnetically ordered states neighboring the AVL fall into two 
categories: conventional XY spin-ordered phases and "super- 
solids," which additionally develop order. 



1. XY spin-ordered states 

Consider the addition of a mass term mAi^^ which is 
favored by a large negative U2 in Eq. ( I63> . Microscopically, 
this mass can be identified with a staggered vortex chemical 
potential that causes the vortices to preferentially occupy one 
of the two sublattices of the honeycomb. The resulting state is 
the vortex "charge density wave" (CDW) shown in Fig.|8ja), 




FIG. 8: Vortex charge density waves (CDWs) proximate to the AVL 
in the isotropic limit J — J' , along with the corresponding spin 
structures. Vortices preferentially occupy filled honeycomb sites. On 
the right side the satisfied bonds of the triangular lattice are solid, 
dashed lines represent less satisfied bonds, and filled circles denote 
sites whose spins fluctuate around zero mean. 



where the vortex density is enhanced on the filled sites and 
depleted on the open sites. 

To identify the corresponding spin structure, first recall that 
the leading spin correlations in the AVL occur at wave 
vectors ±Q and Ki.2,3 of Fig.0 It is natural, then, to expect 
that magnetically ordered descendants of the AVL will involve 
condensation of 5*+ at these wave vectors. We will assume 
this is the case, and search for the symmetry-equivalent of 
^3 ^ ^3 by considering bilinears involving 5* at these wave 
vectors. The answer is unique (up to an overall sign), and we 
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identify 



(66) 



Since (^1^3x^3) 7^ in this vortex CDW, it follows that the 
spin order can be obtained from (5q) 7^ 0, (S'^q) = (or 
vice versa, depending on the sign of the mass m). This is the 
well known \/3 x V3 spiral state depicted in Fig. EJa). We 
note that our identification in Eq. j66> holds in the anisotropic 
case as well; in this regime an incommensurate spiral results. 

As another example, assume the U4 interaction is strong 
enough that a mass term J^j "^j^j is generated. This mass 
similarly corresponds to a modulated vortex chemical poten- 
tial which drives CDW ordering. With only quartic fermion 
terms, there is a large degeneracy of possible states due to 
the arbitrariness in the relative values of toi,2,3- This de- 
generacy is broken, however, by higher-order terms in the 
action, which select either (I) nrii ^ 0, rrij^i = 0, or (II) 
|mi| = \m2\ = jmsl ^ 0. The resulting vortex CDW's are 
shown in Figs.lSJb) and (c), respectively. 

The spin order in these states can be determined using the 
same logic as above. Here, we identify 



(67) 



Consider case (I) first, with say (/Ci,2) = and (/C3) 7^ 0. 
Equation (I67> then implies that {S^,^) = and (5*^^ ^) ^ 
g«¥'i,2^ which yields (/C3) ~ sm{ipi — Lp2)- To maximize 
K/Cs)!, which is energetically favored by the large U4 inter- 
action, the phases ipi and ip2 are chosen to differ by tt/2. 
The resulting spin order is shown on the right-hand-side of 
Fig- Etb)- Now let's consider case (II), where |(/Ci)| = 
|(/C2)| = |(/C3)| ^ 0. Here we take {S+^) - e^'^^ yielding 
(/Ci) ~ sin((/32 — 'Ps), etc. In this case the phases 1^1,2.3 must 
differ by either tt/3 or 27r/3, which leads to the spin order 
illustrated in Fig. ^c). On the right-hand-side of Fig. |SJ the 
solid lines indicate satisfied bonds, while the filled circles de- 
note triangular lattice sites with spins fluctuating around zero 
mean. 

A similar analysis can be used to identify the states arising 
from spontaneously generated mass involving the enhanced 
bilinears Vj. Such mass terms give rise to modulated nearest- 
neighbor hopping amplitudes for the vortices, and drive vortex 
"valence bond solid" (VBS) order Since Vj carries momen- 
tum Kj + Q, the spin structures corresponding to these VBS 
phases involve condensation of both 5q and . At present 
it is unclear what order is driven by IC'j mass terms, due to 
the fact that they arise from second-neighbor vortex hopping, 
which does not have a clear interpretation for the spin system. 



2. Supersolids 

Finally, consider a mass term m[e^'^ M^gg + H.c] generated 
by a large M3 interaction. Microscopically, this mass induces 
both modulations in the nearest-neighbor vortex hopping am- 
plitudes and modulations in the gauge flux piercing the hon- 
eycomb plaquettes. The degeneracy in the phase 7 is lifted by 




FIG. 9: Supersolids neighboring the AVL. On the left side, the di- 
rection of induced gauge flux is specified by a ± sign, and the hop- 
ping amplitudes are dominant along bold links of the honeycomb. 
The "solid" ordering pattern for (S*^) follows the pattern of induced 
gauge flux, while the "superfluid" order for {S^) is shown on the 
right side. Filled circles denote sites with {S^) — 0. The overall 
spin structure can be viewed as a spiral state, tilted into the {S^ ,S^) 
plane. 



higher-order terms in the action, which select either 7 = n7r/3 
or 7 = {2n + l)7r/6, where n is an integer. The vortex states 
for these two cases are shown in Figs. |9ja) and (b), respec- 
tively. In the figure, the direction of induced flux through a 
given plaquette is indicated with a ± sign, and the bold hon- 
eycomb links have the dominant hopping amplitudes. In Fig. 
|9ja), the induced flux is twice as large on "+" plaquettes; in 
(b), the induced flux is equal and opposite on "±" plaquettes 
and vanishes on others. 

These lattice-scale gauge flux modulations signify the onset 
of ordering in the spin system, 



(5") -cos(Q-r + 7). 



(68) 



Since there is a gapless photon in these states, the continu- 
ous U(l) spin symmetry is also broken. Hence the in-plane 
spin components order as well, so that these states are ex- 
amples of "supersolids". Using the symmetry of the vortex 
phases and the uncertainty principle as a guide, the simplest 
assumption for the 5*+ order is shown on the right side of Figs. 
|9ja) and (b). The filled circles in the figure denote sites with 
(5+) = 0. The in-plane and out-of-plane spin structure can 
be collectively viewed as a coplanar spiral state rotated into, 
for instance, the {S"^ , S^) plane. Both spin patterns exhibit a 
\/3 X \/3 periodicity. The difference is that in Fig.|9ja) spins 
on one sublattice point along the hard axis, while in Fig. 
|9[b) spins on one sublattice point along the 5*^ axis. 
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VI. DISCUSSION 

In this paper we have performed a detailed characteriza- 
tion of a novel "critical" spin liquid, the algebraic vortex liq- 
uid, which arises rather naturally out of a reformulation of 
the easy-plane spin model in terms of fermionized vortex de- 
grees of freedom. Among the most striking predictions for 
the AVL is the behavior of the dynamical spin correlations. 
As a consequence of an emergent global SU(4) symmetry, the 
in-plane spin structure factor exhibits enhanced power 
law correlations with identical exponents at the five inequiv- 
alent momenta shown in Fig. Q Due to the easy-plane char- 
acter of the AVL, the out-of-plane spin correlations behave 
quite differently. The S^^ structure factor exhibits enhanced 
power law correlations only at the spiral ordering wave vec- 
tors (±Q in Fig. 0, and is generally expected to be weaker 
than 5*^ . These nontrivial features in the spin structure fac- 
tor distinguish the AVL from other known spin liquids, and 
should serve as useful characterizations for identifying this 
phase experimentally. 

Our study was partly motivated by the spin- 1/2 triangular 
antiferromagnet CS2CUCI4, whose spin dynamics have been 
explored with neutron scattering.'*'^ Although this material 
develops long range spiral order at low temperatures T < 
Tm — 0.62K, the dynamical structure factor exhibits anoma- 
lous power laws at intermediate energies, both in the ordered 
phase and in a range of temperatures above Tjv. Such power 
law behavior is reminiscent of spin-liquid physics, and several 
scenarios for its origin have been proposed. These include 
physics dominated by one-dimensional chains,'*^'*'* a two- 
dimensional algebraic spin liquid,'^ the Z2 spin liquid,"*^"*^ 
and a quantum critical point scenario,^^ as well as more con- 
ventional explanations such as nonlinear spin wave theory.'*'* 

Possible application of the AVL to CS2CUCI4 has been dis- 
cussed in some detail in Ref. 26. The most intriguing obser- 
vation here is that the dynamical structure factor was found 
experimentally to decay with the same power law near wave 
vectors Ki 2 and Q in Fig.E*" consistent with our expecta- 
tions for the AVL. It is important to keep in mind, however, 
that the dual vortex formulation employed here requires an 
easy-plane U(l) spin symmetry so that vortices exist as sta- 
ble topological excitations. This imparts the AVL with a dis- 
tinct easy-plane character, unlike other theoretical proposals 
which retain full SU(2) spin symmetry. Although there is a 
microscopic easy-plane spin anisotropy in CS2CUCI4 due to 
a Dzyaloshinskii-Moriya interaction, this coupling is fairly 
weak, and it is therefore not clear whether the AVL described 
here can be applied directly (some scenarios are considered 
m Ref. '2?). An interesting possibility is that the AVL has an 
SU(2)-invariant relative which may be relevant for CS2CUCI4, 
though we do not know how to access such a state theoreti- 
cally. One speculation in this respect is that there may exist 
a slave-particle description of the AVL. In particular, the di- 
rect slave-fermion approach often yields critical states dubbed 
algebraic spin liquids (ASLs); such states on the triangular lat- 
tice were explored in Ref. 13. The algebraic vortex liquid is 
not formulated using spinon fields, but predicts critical power 
law spin correlations reminiscent of those in ASLs. It should 



also be said that spinons are strongly interacting in the ASL 
and cannot be thought of as free fields in any sense, and the 
same is true about vortices in the AVL. Unfortunately, so far 
we have been unable to find a spin liquid state on the trian- 
gular lattice that would reproduce all the dominant wave vec- 
tors present in the AVL phase, but such a connection between 
these very different theoretical perspectives remains a tanta- 
lizing possibility. 
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APPENDIX A: TRANSFORMATION PROPERTIES OF THE 
NEGATIVE-ENERGY DIRAC SEA 

This appendix is devoted to obtaining the transformation 
properties of \DS,q){Q\, which are needed for determining 
the quantum numbers carried by the monopole operators in 
QED3. Here \DS,q) is the filled negative-energy Dirac sea 
with a 271(7 Aux insertion, where q — ±1, and |0) is the ground 
state with no added flux. We attack the problem in two stages. 
First, we constrain the transformation properties as much as 
possible using various general relations among symmetries. 
The ambiguities that still remain here are then fixed using nu- 
merical studies of monopole insertions. 



a. General arguments 

Fennionic time reversal and Particle-hole symmetry. By 
examining Table HIl we see that the flux changes sign under 
both Tfcim and C. Hence Tfcrm transforms the filled negative- 
energy Dirac sea with q = +\ into the negative-energy Dirac 
sea with q ~ —1, while C additionally fills the four zero- 
modes since \DS,q) is not half-filled. The ground state |0), 
on the other hand, is an eigenstate of both symmetries. Using 
•^orm = = 1, we can then define the phases of \DS,q) 
such that 

Tfe™ : \DS,q){Q\^\DS,~q){Q\ (Al) 
C : \DS,q){Q\^\{f\^^^^\DS,^q){Q\. (A2) 

Aa 

Both \DS, q) and |0) are expected to be eigenstates under 
the remaining symmetries in Table IVl all of which leave the 
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flux invariant. Quite generally, we then have 

5'^"\DS, q){0\ 



R-JT 



\DS,q){0\ 
\DS,q){0\ 
\DS,q){0\ 
\DS,q){0\ 



■-\DS,q){0\ 
^'\DS,q){0\, 



(A3) 



where the last line holds only in the isotropic limit. We will 
now examine the general constraints on the above eigenvalues. 

Inversion. First, one can show that the phases 0^ 
must be independent of q by using the commutation re- 
lation [Rtt , Tform] = when acting on half-filled states, 
which are gauge invariant. For example, it follows from 
[R^,Tf,,^]F^ \DS,q) = that 09 = 0-^. Furthermore, 



R^ = 1 implies that e 



lei — 



Ctt = ±1, so we have 



R^ : \DS,q){0\^C.\DS,q){0\. 



(A4) 



The sign is in principle fixed, but can not be determined 
using general relations alone. 

Modified reflection. Similarly, it follows from the commu- 
tation relations [JZx , C] = [TZx , Tform] = (on physical states) 
that 



ilx : \DS,q){0\^iCx\DS,q){0\. 



(A5) 



The sign Cx = ±1 is also fixed, but can not be determined 
from this general analysis. 

Translations. We constrain the phases Of ^ by first as- 
suming the following operator relations hold when acting on 
gauge-invariant states. 



T2R-K — RirT2 ^, 

T2'Tix — Ti-xTiT2 , 



(A6) 
(A7) 



since the left and right sides transform the lattice identically. 
The first relation implies e*^2 = ±1, while using the second 
we conclude e*^i = —1. We can fix the former sign by now 
specializing to the isotropic limit. Here, we have an additional 
symmetry relation. 



TiRtt/S — Rtt/3T2 



(A8) 



that holds when acting on physical states. From this we obtain 
_ git?2 — _2 in the isotropic limit. By continuity, we 
assume this carries over in the anisotropic limit as well so that 



T,^2:\DS,q){0\^-\DS,q){0\. 



(A9) 



Rotations (isotropic limit). The commutation relation 
[i?7r/3 , Tform] = on physical states implies that e 
g~'^7r/3 Moreover, the relation = R-n together with 



x/3 — 



commutation with particle-hole symmetry yields e '/^ 
CttC'^*''^/"^. Thus, we have 



i?^/3 : \DS,q){0\ Cxe2'9-/3|i5^,5)(0|. 



(AlO) 



We have now arrived at the transformation properties listed 
in Eq. j59l l. As discussed in Sec. llV C^ determining the am- 
biguities in (^7r and C,x that arose above is crucial for under- 
standing the in-plane spin correlations in the AVL. We will 
attempt to sort out these uncertainties by appealing to numer- 
ics, discussed below. 



b. Numerical diagonalization 

For convenience, we specialize to the isotropic limit J — J' 
for the remainder of this appendix. Consider the mean-field 
Hamiltonian M6\ generalized to include arbitrary flux inser- 
tions. 



-t 



(xx'> 



H.c.]. (All) 



Here is the gauge field giving rise to a background tt flux as 
before, while 5a gives rise to any added flux. We have numer- 
ically diagonalized the above Hamiltonian on finite systems 
to obtain the spectrum of single-particle energies and wave 
functions. Specifically, we considered lattices composed of 
Airings concentric "rings" of honeycomb sites. For example, 
Figs.|8land|9]illustrate systems with Airings = 2 and 3, respec- 
tively. System sizes up to A'rings = 17, consisting of 1734 
lattice sites, were studied. Flux insertions varying from to 
27r were taken to be uniformly spread within the first several 
innermost rings. A variety of boundary conditions were used, 
namely, open boundary conditions; "Klein bottle" boundary 
conditions, where one connects boundary sites with coordina- 
tion number 2 and their inversion counterparts; and modifica- 
tions of the latter, where one connects only a subset of such 
boundary sites and their inversions. (Note that open boundary 
conditions are problematic when A'lmgs is odd, because with 
no added flux there are two degenerate zero-energy modes and 
therefore no unique ground state |0). Other boundary condi- 
tions mentioned above give a unique ground state as desired.) 
Such geometries break translational invariance, but are par- 
ticularly convenient for addressing the symmetries of interest 
here. 

With the single-particle wave functions in hand, one can 
explicitly construct the states |0) and \DS, +1). The ground 
state |0) is simply built from all negative-energy wave func- 
tions. More care must be taken, however, in constructing 
+1). In a finite system, the four quasi-localized "zero- 
modes" — which are excluded from this state — are split away 
from zero energy, two above zero and two below. Identifying 
these quasi-localized modes is essential, particularly since the 
boundaries can introduce additional "edge" modes near zero 
energy. These modes can be distinguished by the behavior of 
their wave functions. Most of the probability weight lies near 
the flux insertion for the quasi-localized modes, whereas the 
dominant weight for the edge modes occurs near the bound- 
ary. A useful diagnostic for this comparison is the "ring partic- 
ipation" Pn(0), which for a particular wave function </) gives 
the probability weight summed over honeycomb ring n, nor- 
malized by the number of sites in the ring. More explicitly. 
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X 10"^ 




FIG. 10: Ring participation Pn versus the ring index n for the first 
several wave functions above zero energy. The data correspond to 
a system with Airings ~ 14, open boundary conditions, and a 2n 
flux insertion spread within the first four innermost rings. The quasi- 
localized modes (solid curves), whose wave functions are peaked 
near the flux insertion, are clearly distinguishable from other low 
energy modes (dashed lines). 



where A'sitosl"-) is the total number of sites in ring n and i 
is summed over all sites in the ring. Figure [TO] displays the 
ring participation for the first several wave functions above 
zero energy in a system with A'rings = 14, open bound- 
ary conditions, and 2tt flux inserted within the first four in- 
nermost rings. This illustrates the clear difference between 
the quasi-localized modes (solid lines) and other low-energy 
states (dashed lines). In most cases observed this distinction 
allows one to identify the former, which are the modes of in- 
terest. Once these have been located, the state \DS, +1) can 
be built out of the remaining negative-energy wave functions. 
The inversion eigenvalues for |0) and \DS,+1) are then 



simply given by the product of inversion eigenvalues for the 
single-particle wave functions contributing to these states. We 
find that the parity under inversion for |0) and \DS, +1) de- 
pends on both the system size and boundary conditions. How- 
ever, in all cases where the quasi-localized modes can be 
clearly resolved, the product of the inversion eigenvalues of 
|0) and \DS, +1) — which gives the sign (fi — is geometry in- 
dependent. In particular, we find (ij ~ —1. 

Some insight into this result can be obtained by viewing 
the 27r flux as being inserted adiabatically. Numerically, this 
is achieved by ramping the added flux from to 2tt in sev- 
eral stages and monitoring the evolution of the energy levels 
during the insertion. In all cases studied, no levels cross zero 
energy during the evolution (though we do not have an argu- 
ment for why this is the case). This implies that the quantum 
numbers for the state with all negative-energy modes occupied 
are unchanged by the 2tt flux insertion. Moreover, we ob- 
serve that the two quasi-localized modes with energy slightly 
below zero always have opposite inversion eigenvalues. Con- 
sequently, the states |0) and \ DS, +1) must also have opposite 
inversion parity. Why the quasi-localized modes split in this 
way is unclear at the moment, but would be useful to explore. 

Obtaining the sign from numerics is more subtle due 
to the antiunitarity of TZx- A more useful symmetry to ex- 
amine is the unitary operation TZ'^ = TZxCTfcmi, which has 
eigenvalues ±1. We will use the eigenvalues of TZ'^ to back 
out the sign (x- The fact that we find no zero-energy level 
crossings provides a useful shortcut to this end (but is not 
necessary). Again, the quantum numbers of the state with 
all negative-energy modes filled are then conserved under a 
27r flux insertion. In particular, the ground state |0) and 
the q — 1 Dirac sea with the two negative-energy quasi- 
localized modes filled must have identical eigenvalues un- 
der both TZ'x and R^/^. There are just two candidates for 

the latter state, since only (i^ +i ^ -^"1 +1) and 
(fI^^ + fI^^ + fI^^)\DS,+1) have the same rotation 
eigenvalue as the ground state |0). Furthermore, under TZ'^ 
both candidates have eigenvalues that differ from |0) by — f^^;. 
Hence in either case we conclude that Q = ^ 1 ■ 
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TABLE VII: Transformation properties of the 48 bilinears which are components of conserved currents at the AVL fixed point and thus do not 
have enhanced correlations. In the first column cr± = cr"" ± ia'" . 







Tlx 




c 




Rt^/3 (isotropic limit) 




P 


P 


P 


-p 


P 


P 


ip^ a±'tl) = G± —> 


G± 


G± 


-G± 


-G± 


-Gt 


e±^tG± 




M"" 


M^" 




M"" 








Mi" 


Ml' 




Mi" 


M^" 


-e±**Mi" 


^/)f(T" +iA("r'')V' = N ^ 




N 


Art 


N 


iV+ 


A^t 






N± 


-Art 


N± 


-iVi 


e*'*Ar^ 


■0 V-r^V' = iCi^ 




K.2 


/Ci 


/Ci 




/C3 








/C2 


IC2 


-/C2 


/Ci 


V'V^V' = ^3 — » 






JC3 


JC3 


-/C3 


/C2 


lX'T''cr±%l) = iCi± — » 






— /Ci± 


lCi± 


/Clip 




ip^ fi^T'^a±ip = ^2± — > 




lCi± 


— /C2± 


K,2± 


A^2^ 




V'V^O'iV' = ^3± — » 




K,3± 


-^3± 


^3± 


^3t 


e±^t/C2± 






—K,2 


-)C[ 








^/-V'-r"^ = /c^i ^ 


6-^^,2 




^IC'2 


~1C'2 


IC'2 










-K', 








t/'>"r''a±(,' = IC'i- 
















e K,.2± 


—K.'i± 


^2± 


-^2± 


— ^2=F 






<-' /C-.ii 












(/' (//'■ + = Pi 




V2 




-•Pi 


-•P| 


PJ 






Vi 


vl 


-P2 




v\ 






V3 


n 


-P3 




v\ 


+ ifi''T')a±^ = Pi± 




V2± 




-Pl± 


vL 








Vi± 


-vL 


-V2± 


vL 








V3± 


-vL 


-V3± 


vL 


e±'*P2% 



L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. Sengupta, 
Phys. Rev. B 71, 144508 (2005). 

L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. Sengupta, 
Phys. Rev. B 71, 144509 (2005). 

A. A. Burkov and L. Balents, Phys. Rev. B 72, 134502 (2005). 
J. Alicea, O. I. Motrunich, M. Hermele, and M. P. A. Fisher, Phys. 
Rev. B 72, 064407 (2005). 
^ J. Alicea, O. I. Motininich, and M. P. A. Fisher, Phys. Rev. Lett. 
95, 247203 (2005). 

W. Zheng, J. O. Fjaerestad, R. R. P Singh, R. H. McKenzie, and 

R. Coldea, cond-mat/0506400 (unpublished). 

R. G. Melko, A. Paramekanti, A. A. Burkov, A. Vishwanath, D. N. 

Sheng, and L. Balents, Phys. Rev. Lett. 95, 127207 (2005). 

D. Heidarian and K. Damle, Phys. Rev. Lett. 95, 127206 (2005). 
'° S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 (2005). 
" M. Boninsegni and N. Prokof'ev, Phys. Rev. Lett. 95, 237204 

(2005). 

W. Zheng, R. H. McKenzie, and R. R. P Singh, Phys. Rev. B 59, 
14367 (1999). 

T. W. Appelquist, M. Bowick, D. Karabali, and L. C. R. Wijew- 

ardhana, Phys. Rev. D 33, 3704 (1986). 

W. Rantner and X.-G. Wen, Phys. Rev. B 66, 144501 (2002). 

M. Franz, Z. Tesanovic, and O. Vafek, Phys. Rev. B 66, 054535 

(2002). 



M. Franz, T. Pereg-Bamea, D. E. Sheehy, and Z. Tesanovic, Phys. 
Rev. B 68, 024508 (2003). 
" K. Kaveh and I. R Herbut, Phys. Rev. B 71, 184519 (2005). 
M. Hermele, T. Senthil, and M. P. A. Fisher, Phys. Rev. B 72, 
104404 (2005). 

S. J. Hands, J. B. Kogut, and C. G. Strouthos, Nucl. Phys. B 645, 
321 (2002). 

S. J. Hands, J. B. Kogut, L. Scorzato, and C. G. Sttoutiios, Phys. 
Rev. B 70, 104501 (2004). 
R. Jackiw, Phys. Rev. D 29, 2375 (1984). 
'"'^ V. Borokhov, A. Kapustin, and X. Wu, J. High Energy Phys. 11, 
049 (2002). 

"■^ M. Bocquet, F. H. L. Essler, A. M. Tsvelik, and A. O. Gogohn, 
Phys. Rev. B 64, 094425 (2001). 



47 



O. A. Starykh and L. Balents (unpublished.). 



44 

C. H. Chung, J. B. Marston, and R. H. McKenzie, J. Phys.: Con- 
dens. Matter 13, 5159 (2001). 

C. H. Chung, K. Voelker, and Y. B. Kim, Phys. Rev. B 68, 094412 
(2003). 



S. V. Isakov, T. Sentiiil, and Y. B. Kim, cond-mat/0503241 (un- 
published). 

M. Y. Veillette, A. J. A. James, and F. H. L. Essler, cond- 
mat/0506667 (unpublished). 



